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ABSTRACT 

The y-ray pulsar PSR B 1706-44 and the adjacent supernova remnant (SNR) candidate G 343. 1-2.3 were observed by H.E.S.S. during a dedicated 
observation campaign in 2007. As a result of this observation campaign, a new source of very-high-energy (VHE; E > 100 GeV) y-ray emission, 
| HESS J1708-443, was detected with a statistical significance of 7 cr, although no significant point-like emission was detected at the position of 
the energetic pulsar itself. In this paper, the morphological and spectral analyses of the newly-discovered TeV source are presented. The centroid 
of HESS J 1708-443 is considerably offset from the pulsar and located near the apparent center of the SNR, at Otaooo = 17 h 08 m ll s ± 17 s and 
<5j20oo = -44°20' ± 4'. The source is found to be significantly more extended than the H.E.S.S. point spread function (~0. 1°), with an intrinsic 
Gaussian width of 0.29° ± 0.04°. Its integral flux between 1 and 10 TeV is ~ 3.8 x 10~ 12 ph cnT 2 s -1 , equivalent to 17% of the Crab Nebula flux 
in the same energy range. The measured energy spectrum is well-fit by a power law with a relatively hard photon index T = 2.0 ± 0. l Btat ± 0.2 sys . 
Additional multi-wavelength data, including 330 MHz VLA observations, were used to investigate the VHE y-ray source's possible associations 
with the pulsar wind nebula of PSR B 1706-44 and/or with the complex radio structure of the partial shell-type SNR G 343. 1-2.3. 
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1. Introduction 

The energetic pulsar PSRB1706-44 (also PSR J1709-4429) is 
one of the first pulsars from which pulsed emission was de- 
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tected not only in the radio dJohnston et al.|[l992h and in X-rays 
dGotthelf et al.ll2002l). but also i n high-energy (HE; E ~ GeV) 
y-rays dSwanenburg et al.lll981l) . Th e pulsar was first dete cted 
in a high-frequency radio survey by [Johnston et al.l (1 19921) and 
has a spin period P — 102 ms, a characteristic age r c = 17 500 
yr, and a spin-down luminosity E = 3.4 x 10 36 ergs~'. It belongs 
to the class of relatively young and powerful pulsars, of which 
the Vela Pulsar is the most prominent example in the south- 
ern hemisphere. The putative wind nebulae of these pulsars are 
prime candidates for being sources of very-high-energy (VHE; 
E > 100 GeV) y-rays. A bright, HE y- ray source, 2CG 342-02 
was discovered by the COS-B satellite dSwanenburg et al.ll 198 lb 
and later found to be positionally coincident with the radio pul- 
sar. The y-ray source was firmly associated with PSRB 1706-44 
after EGRET (the Energetic Gamma Ray Experiment Telescope, 
onboard the Compton Gamma-Ray Observatory) observed pul- 
sations from 3EG J 17 10-4439 (also EGRJ17 10-4435) which 
match ed the period seen in the radio waveband (iThompson et alj 
Il992l) . More recently, the pulsar has been detected at GeV en- 
ergies by the latest generation of spaceborne HE y-ray detec- 
tors: by AGILE (Astro rivelatore Gamma ad Immagini LEggero) 
as 1AGL J 1709-4428 dPittori et al.ll2009l) and by Ihe Fermi/LAT 
(Larg e Area Telescope) as 1 FGL J 1 709 .7-4429 dAbdo et alj 
1201(1 . 

Radio observations of PSRB1706-44 reveal the presence 
of a synchrotr on nebula, with an extension of -3', surround- 
ing the pulsar dFrail et al.lll994t iGiacani et aT]|2001l) . The ob- 
served polarization and the flat spectrum, with a flux density 
spectral slope a = 0.3 (where the flux density S oc v~ a ), 
suggest it is a pulsar wind nebula (PWN). However, the im- 
plied conversion efficiency from spin-down energy to radio flux 
of -2 x 10~ 6 would be the lowest of any known radio PWN 
dGiacani et alj 1200 lb . Observations by the X-ray telescopes 
onboard ROSAT (Roentgen Satellite) and ASCA (Advanced 
Satellite for Cosmology and Astrophysics) reve al that the neb- 
ula is also visible in X-rays dFinlev et al.lll998l) . The morphol- 
ogy of the PWN was mapped in detail at arcsecond scales uti- 
lizing the superior re solution of the Chandra X-ray Observatory 
dRomani et al.ll2005l) . The X-ray analyses suggest the presence 
of a diffuse X-ray PWN, with a spectral index of 1 .77, which sur- 
rounds a more complex structure consisting of a torus and inner 
and outer jets. The diffuse X-ray PWN has a radius of 1.8' and 
also exhibits a fainter, longer extension to the West. The pres- 
ence of non-deformed X-ray jets is consistent with the pulsar's 
low apparent s peed, v — 89 km s 1 , a s deduced from scintillation 
measurements dJohnston et al.lll998l) . 

The pulsar PSRB 1706-44 is also located at the southeast 
end o f an incomplete arc of radio emission dMcAdam et al.l 
119931) . which has been suggested to be the partial shell of a 
faint supernova remnant (SNR G 343. 1-2.3). The arc is embed- 
ded in weak diffuse radio emission, which is present both in- 
side and outside of the arc ([Frail et al 1 11994 . Polarization mea- 
surements suggest that this diffuse emission i s associated with 
synch rotron radiation from the SNR itself dDodson & Golapl 
120021) . T he SNR has not been detected in any othe r waveband 
(see e.g. iBecker et al.Hl995b lAharonian et alj|2005al) . There are 
various estimates of the dis t ance to the pulsar, ran ging from 
1.8 koc dJohnston et al I[l992l: iTavlor & Cordesll 19931) to 3.2 kpc 
dKoribalski et al.l 1995b . The distance 2.3 + 0.3 kpc, derived from 
the dispersion meas ure and the most recent Galactic free electron 
distribution model dCordes & Lazioll2002l) . is adopted through- 
out this paper. This distance is compat ible with the less relia ble 
2 - D distance of -3 kpc for the SNR dMcAdam et al.lll993l) . 



The possible physical association between PSRB 1706-44 
and G 343. 1-2.3 has been questioned based on the differing 
age and distance estimates for the SNR and pulsar (see Sec. 
4.2 and 4.3, respectively) and the lack of visible interaction. 
Furthermore, if the pulsar originated at the apparent center of the 
SNR, then its inferred velocity (-700 k m s" 1 ) is incompatible 
with it s scintillation velocity (89 km s -1 ). lBock & Gvaramadzel 
d2002l) suggested a scenario where an off-center cavity explosion 
could relax the restrictions on the inferred velo city and invali- 
date t he age estimate for the SNR of -5 000 yr dMcAdam et al.l 
[19931) . which is based on a Sedov-Taylor model. In this sce- 
nario, PSRB 1706-44 and G 343. 1-2.3 are physically associ- 
ated; however, the radio arc is not identified with the SNR shell, 
but rather with the former boundary of the wind-blown cav- 
ity that was overtaken a nd compressed by the expanding SNR 
dDodson & Golapll2002l) . 

In the VHE domain, both the pulsar and the SNR have been 
observed using ground-based, imaging atmospheric-Cherenkov 
telescopes (IACTs). The findings of the various observations are, 
however, not fully consistent (see Table[TJ. The CANGAROO 
(Collaboration of Australia and Nippon (Japan) for a Gamma 
Ray Observatory in the Outback) Collaboration reported the 
detection of steady emission, coincident with the pulsar, using 
the 3 8 m CANGAROO-I telescope in 1992-1993 kifune et all 
Il995l) . They measured an integral flux F(> 1 TeV) « 1 x 
10~ n ph crrT 2 s" 1 , equivalent to -44% of the Crab Nebula fluxQ 
However, the CANGAROO Collaboration recently undertook a 
comprehensive re-analysis of their archival CANGAROO-I data 
and no longer find a signal; instead, they calculate an upper 
limit (UL; here, at 95% confidence level and assuming a spec- 
tral index of -2.5) to the integral flux of F(> 3.2 TeV) < 8.0 x 
10" 13 ph cirr 2 s 1 (-24% Crab) dYoshikoshi et al.ll2009l The 
4-m BIGRAT (BIce ntinnial Gamma RAy Telescope) telescope 
dRowell et al.l 1 19981) also observed the pulsar in 1993-1994 
and reported a compatible UL. Observations in 1 996 with the 
Durham Mark 6 telescope dChadwick et al.l [l 9981) appeared to 
confirm the earlier CANGAROO-I detection, with a reported 
integral flux that was compatible within the large systematic 
uncertainties (+30% for CANGAROO-I and +50% for the 
Mark 6). Further observations with the CANGAROO-II tele- 
scopes in 2000-2001 a gain seemed to validate the detection 
dKushida & et al.1 120031) . However, when the H.E.S.S. (High 
Energy Stereoscopic System) Collaboration observed the pul- 
sar in 2003 during its commissioning phase, they did not de- 
tect any significant VHE y-ray emission from PSRB 1706-44 
or its vicinity. The derived UL (99% confidence level) on the 
integral flux from an extended region encompassing the SNR 
was fo und to be F(> 0.5 TeV) < 3.5 x 10~ 12 ph cirr 2 s _1 (-5% 
Crab) dAharonian et al 1 l2005al) . in stark disagreement with all 
of the previous findings (see also the Appendix). Shortly there- 
after, preliminary analysis of stereo observations with the 4x10- 
m CANGAROO-III telescope array also disagreed with the ini- 
tial CANGAROO-I detection and resulted in an UL at the pul- 
sar position (95% confidence l evel) of F(> 0.6 TeV) < 5 x 
10~ 12 ph cirr 2 s 1 (-10% Crab) dTanimori & et al.ll2005l) . which 
agreed with the H.E.S.S. results at the time. 

In 2007, H.E.S.S. followed up on their initial result by carry- 
ing out additional dedicated observations of the pulsar, now uti- 
lizing the superior sensitivity of the fully-operational H.E.S.S. 
telescope array. This campaign resulted in the discovery of 
extended emission from the vicinity of PSRB 1706-44 and 



1 The Crab Nebula spectrum published in Aharonian et al. (2 006af) is 
adopted as the reference Crab spectrum throughout this paper. 
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Appendix A 
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PSRB 1706-44 
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F(> 1 TeV) = (3.0 ± 0.6) x 10~ 12 


Enomoto et al. (2009) 
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CANGAROO-III 
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Enomoto et al. (2009) 
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H.E.S.S. 


PSRB 1706-44 
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Sect. 3 


2007 


H.E.S.S. 


G 343. 1-2.3 center 
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Table 1. Summary of the efforts to observe PSR B 1706-44 in the VHE y-ray domain. The CANGAROO upper limits (ULs) are at 
a 95% confidence level (CL), the BIGR AT UL is at 3 cr ( 99.7% CL), and the H.E.S.S. ULs are at a 99% CL. The CANGAROO-I 
integral flux based on the 1992 dataset (lOgio et al.|[l9 93) likely suffered from systematics similar to those that affected the 1993 
data, which has since been re-analysed along wit h an analysis of the pre viously unreleased 1994-1998 CANGAROOT data; the 
ULs assume a Crab-like spectral index T = -2.5 (lYoshikoshi et al]l2009l) . Only the latest, revised results are shown in this table ; 
see Sect. 1 for further discussion. The BIGRAT UL is subject to an additional + 50% systematic uncertainty dRowell et al.lll998h . 
The integral flux from t he 2000-2001 CANGAROO-II data analysis was not disclosed but was claimed to confirm previous results 
dKushida&etalJl2003h . The 2003 H.E.S.S. ULs are based on the re-analysis presented in Appendix A; the first UL (row 10) 
assumes T = -2.5 for comparison to the CANGAROO ULs, while the second UL (row 11) assum es T - -2.0 for compa rison to the 
2007 H.E.S.S. detection. The CANGAROO-III fluxes are from the ON-OFF analysis presented in lEnomoto et al.l (120091) . The 2007 
H.E.S.S. results are described in Sect. 3, where the center of G 343.1-2.3 is also defined; the point-source UL assumes F = -2.5. 



G 343 .1-2.3, with preliminary results published in iHoppe et al.l 
(120091) . The latest results from CANGAROO-III also indicate 
the presence of an extended source of VHE y-ray emission from 
the vicinity of the pulsar, although the spectrum and morphol- 
ogy of the emission vary consider ably depending on th e method 
used for background subtraction (lEnomoto et al.ll2009l) . For ex- 
ample, integrating within 1 .0° of the pulsar position and using 
an ON-OFF background method (see Sect. 2.3), they find a Crab 
Nebula-level integral flux. In this paper, we present new VHE 
data on PSRB 1706-44 and G 343.1-2.3 which was obtained 
during H.E.S.S.'s 2007 observational campaign. 



2. H.E.S.S. Observations and Analysis 

2.1. The H.E.S.S. Telescope Array 

H.E.S.S. is an array of four IACTs, dedicated to the observation 
of VHE y-rays. The array has been operating since December 
2003 in the Khomas Highlands of Namibia; its location in the 
southern hemisphere (23° 16' 17" S) allows observations of the 
inner Galaxy at reasonably low zenith angles. Each telescope is 
equipped with a tessellated, spherical mirror with a total area of 
107 m 2 and a camera comprised of 960 photomultiplier tubes, 
covering a field-of-view (FoV) 5° in diameter. The telescopes 
are situated on a squar e with sides of 120 m length and operated 
in stereo trigger mode dFunk et al.ll2004l) . which requires at least 
two telescopes to trigger the detection of an extended air shower 
(EAS). This stereoscopic approach results in an angular reso- 
lution < 0.1° per event, an energy resolution of ~15% (on av- 
erage) , and an improved background rejection (lAharonian et alj 
l2006ah . The H.E.S.S. array can detect point sources near zenith 
at flux levels of ~1% of the Crab Nebula flux with a statisti- 
cal significance of 5 cr in 25 h of observations, o r less if ad- 
vanced techniques are used for EAS image analysis (lAcero et alJ 



120091) . Its large FoV and off-axis sensitivity not only make it 
ideally suited for surveying the Gal actic Plane (lAharonian et al.l 
I2005bl I2006ct IChaves et all 120081) . but also for studying ex- 
tended sources like HESS J1708-443. 



2.2. VHE y-ray Observations 

The region of interest, which includes PSRB 1706-44 and 
SNR G 343.1-2.3, was observed with the full four-telescope 
H.E.S.S. array in 2007. The observations were dedicated to 
search for VHE y-ray emission from the pulsar and were there- 
fore taken in wobble mode centered on its positio n in the ra- 
dio (g ?0 no = 17 h 09 m 42.73 s , c5 2 ooo = -44°29'08.2"; IWang et al.l 
(I2000h ). In this observation mode, the array is pointed toward 
a position offset from the source of interest to allow simul- 
taneous background estimation. Observations of 28-min dura- 
tion were taken, alternating between offsets of ±0.7° in dec- 
lination and right ascen sion. After standard quality selection 
(lAharonian et al.l 1200 6a) to remove data affected by unstable 
weather conditions or hardware-related problems, the total live- 
time of the dataset is ~15 h. The zenith angle of the observations 
ranges from 20° to 30°, with a mean of 24°. We only use data 
from the 2007 observations of PSRB 1706-44, because at that 
time the full four-telescope array was in operation along with 
the central stereo trigger system, resulting in a higher sensitiv- 
ity com pared to earlier observations in 2003 (lAharonian et al.l 
l2005al) when H.E.S.S. was in its commissioning phase, with only 
two telescopes and no central trigger (see also the Appendix). 



2.3. Analysis Methods 

The dat aset was ana lyzed using the Hillas second moment 
method (iHillasI 119851) and t he H.E .S.S. standard analysis de- 
scribed in Aharonian et al . (2006a). For y-hadron separation, 
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17 h 10 m 17 h 05 m 17 h 12 m 17 h 10 m 17 h 08 m 17 h 06 m 17 h 04' 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 1. Left: Image of the VHE y-ray excess (in units of y-rays arcmin~ 2 ) from HESS J1708-443, smoothed with a 2D Gaussian 
with a width <x = 0.10°. The blue-to-red color transition is chosen to reduce the appearance of features which are not statistically 
significant. The white cross is located at the best-fit position of the center-of-gravity of the emission and its size represents the 
statistical error of the fit. The small and large dotted white circles, labeled A and C, respectively, denote the regions used for spectral 
analysis. The a priori defined Region B, from which the detection significance was calculated, is represented by a dotted green 
circle. The three regions are summarized in Table |2] The position of the pulsar PSR B 1706-44, at the center of region A, is marked 
by a square. The inset (bottom-left corner) shows the point-spread function of the H.E.S.S. telescope array for this particular dataset, 
smoothed in the same manner as the excess image. Radio contours of constant intensity, as seen at 330 MHz with the Very Large 
Array (VLA), are shown in green. The radio data were smoothed with a Gaussian of width <x = 0.03°. The local maximum in the 
radio contours at the center of the image is largely due to PMN J1708-4419, an extragalactic object seen in projection (see Sect. 
4.3). Right: Gamma-ray excess in quadratic bins of 0.175° width. The upper number in each bin is the excess summed within this 
bin, and the lower number is the corresponding statistical error. The blue contours correspond to a smoothed excess of 0. 14, 0.17, 
and 0.21 y-rays arcmin~ 2 , taken from the image on the left. The red-rimmed bin is centered on the pulsar position. Note the different 
field-of-view used in the two figures. 



hard cuts were used, which require a minimum of 200 pho- 
toelectrons (p.e.) to be recorded per EAS image. Compared to 
standard cuts (80 p.e.), this relatively strict requirement results 
in better background rejection and an improved angular reso- 
lution but also in an increased energy threshold (560 GeV for 
this dataset). The time-dependent optical response of the system 
was estimated from the Cherenkov light of single muons passing 



close to the telescopes ( Bolz 20041). Thre e different background 
estimation procedures dBerge et al. 2007]) were used in this anal- 



ysis. 



For 2D image generation, the ring background method 



dBerge et al.ll2007l) was used with a mean ring radius of 0.85 c 
Since this method includes an energy-averaged model for the 
camera acceptance to account for the different offsets of the sig- 
nal and background regions from the camera center, it was not 
used fo r spectral extractio n. The reflected region background 
method dBerge et al. 2007) was instead used to measure the flux 
from the pulsar position. 

Since the observations of PSRB 1706-44 were performed 
in wobble mode (see Sect. 2.2), half are actually pointed in- 
side the extended emission from HESS J1708-443, which was 
not known to exist at that time. Therefore, for spectral extrac- 
tion from extended regions which also enclose the pointing po- 
sitions of the telescopes, the bac kground was estim ated using 
the ON-OFF background method dBerge et al.ll2007l) . where off- 
source (OFF) data taken is taken from extragalactic regions of 
the sky where no y-ray sources are known. To match the observ- 
ing conditions between on-source (ON) and OFF data, the two 



observations had to be taken within six months of each other 
and at similar zenith angles. The ON-OFF ba ckground method 
was al so used for the analysis of Vela Junior dAharonian et al.l 
l2005cl). The norm alization between ON and OFF observations 
dBerge et alj|2007l) was calculated from the total event number 
in the two observations, excluding regions with significant VHE 
y-ray signal. The background is thus normalized in an approx- 
imately ring-shaped region (depicted in Fig. 3) with inner ra- 
dius 1.0° and outer radius 2.5°, excluding a sm all region which 
overl aps the known source HESS J1702-420 dAharonian et al.l 
120081 With this background normalisation, the analysis is ob- 
viously only sensitive to a localized excess of y-rays but not to 
emission which would be more or less uniform across the entire 
H.E.S.S. FoV. 



3. Results 

Two different circular regions were defined a priori in order 
to reduce the number of trials during a search for statistically- 
significant VHE y-ray emission. Since other IACTs have re- 
ported point-like emission from the pulsar position, one of these 
regions, hereafter Region A, is centered at this position and has 
a radius of 0. 10°, which is the standard radius used to search for 
point sources in the H.E.S.S. standard analysis. The second re- 
gion, hereafter R egion B, is identi c al to th e region referred to as 
the Radio arc in lAharonian et al. I J2005ah : it is centered at the 
approximate apparent center of the radio arc (0:2000 = 17 h 08 m , 
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Region 


Center 

Q '2000 ^2000 


Radius 
[°] 




N off 


a 


Excess 


Significance 
[o-] 


Integral flux (> 0.6 TeV) 
[ph cm" 2 s~'] 


A 


17 h 09 m 42.73 s -44°29'8.2" 


0.10 


71 


717 


0.11 


_ 9 4+«-2 


-1.0 


< 3.3 x 10-" 


B 


17 h 08 m -44°16'48" 


0.60 


3180 


2488 


1.06 


543+^7 


7.0 


= 6.5 x 10- 12 


C 


17 h 08 ra ll s -44°20' 


0.71 


4243 


3425 


1.06 


615+ 90 


6.8 


= 6.9 x 10" 12 



Table 2. Event statistics for Regions A, B, and C. The center and the radius of each circular on-source (ON) region is given in 
columns 2-4. For Region A, the background was extracted from off-source (OFF) regions in the same field-of-view, while for 
Regions B and C, it was estimated from observations of separate OFF regions. Due to the smaller extent of Region A, more OFF 
regions could be used, which resulted in a smaller normalization factor a than for Regions B and C. The number of events in the ON 
and OFF integration regio ns, N on and N ff , respectively, are given in columns 5 and 6. The significance (column 7) was calculated 
following the approach of iLi & Mai (|T983). The integral flux (or UL thereof) for each region is given in column 8. Note that the 
statisti cs presented here were ob tained from a dataset comprised only of observations in 2007, which does not overlap with the one 
used in lAharonian et al 1 (l2005ah . 



HESS J1 708-443: Spectrum 




< -i- 



Energy [TeV] 

Fig. 2. Differential energy spectrum of HESS J 1708 -443, ex- 
tracted from Region C (see Table [2j. The solid line shows the 
result of a power-law fit. The error bars denote 1-cr statistical er- 
rors. The bottom panel shows the residuals of the power-law fit. 
Events with energies between 0.6 and 28 TeV were used in the 
determination of the spectrum, and the minimum significance 
per bin is 1 <x. 



£2000 = -44°16'48"; as defined in lAharonian et al.l (l2005al) ) and 
has a radius of 0.60° in order to enclose the entire radio structure. 

No statistically-significant emission is observed from the 
pulsar position (Region A); therefore, an upper limit of 14.8 ex- 
cess y-ray events is calculate d at a 99% confidence level, follow - 
ing the unified approach of iFeldman & Cousinsl (119981) . From 
Region B, however, a clear signal is detected with 543 excess 
y-rays and a significance of 7.0 <x. The measured signal corre- 
sponds to a flux ~13% that of the Crab Nebula above 0.6 TeV. 
Table [2] summarizes the event statistics for Regions A and B. 

Figure Q] (left) presents an image of the VHE y-ray excess in 
the 2° x 2° region around the source, smoothed with a Gaussian 
of width 0.09° to reduce statistical fluctuations. The smoothing 
radius is chosen to be on the same scale as the H.E.S.S. point- 
spread function (PSF; 68% containment radius ~ 0.1°), so that 
resolvable morphological features are largely maintained. The 
emission clearly extends beyond the PSF, which is depicted in 
the lower left inlay of Fig. [TJ (left). Figure [TJ (right) shows the 
number of excess events within the emission region along with 



their statistical errors, in quadratic bins of 0.175° width, without 
smoothing. This figure demonstrates that the current statistics do 
not permit a detailed study of the source morphology. However, 
the lack of a significant VHE y-ray excess at the position of the 
pulsar is clear in both figures. 

The centroid of the new H.E.S.S. source is determined by 
fitting the unsmoothed y-ray excess image with a radially- 
symmetric Gaussian profile (<f> = <f>o e~ r ~/( 2o "~)) convolved with 
the H.E.S.S. PSF (0.07° for this analysis). The centroid of the 
best fit is at a J20 oo = 17 h 08 m l I s + 17 s and feooo = -44°20' ± 4' 
({ ~ 343.06°, b ~ -2.38°). T he pointing precision of the 
H.E.S.S. telescope array is 20" dGillessen et al.l 120051) . which 
adds an additional systematic uncertainty. The combined errors 
are reflected in the size of the cross in Fig.[TJ(left). Consequently, 
the new VHE y-ray source is designated HESS J1708-443. The 
fit also gives the source's intrinsic Gaussian width <x = 0.29° + 
0.04° stat . 

Spectral analyses were performed for two regions, Region 
A, which was introduced above, and Region C, which is cen- 
tered on the centroid (i.e. best-fit position) and has a radius 
of 0.71° (see Table The size of Region C represents an 
~95% enclosure of the excess, chosen as a compromise be- 
tween an optimal signal-to-noise ratio and independence of 
source morphology. Both regions are indicated by dashed cir- 
cles in Fig. [TJ (left). From Region A, an integral flux limit of 
F(> 0.6 TeV) < 3.3 x 10~ 13 ph c m' 2 s" 1 was derived with a 
99% confidence level according to IFeldman & Cousinsl d 1998b . 
assuming that the underlying y-ray spectrum follows a power 
law with photon index T - 2.5, a n index close to that of the Crab 
Nebula (lAharonian et al 1 l2006ah . This upper limit corresponds 
to ~ 1 % of the flux of the Crab Nebula in the same energy range. 

The energy spectrum of the entire source is extracted from 
Region C. Within the large integration circle, 615 excess y-ray 
events were found, corresponding to a statistical significance 
of 6.8 cr (pre-trials). The differential spectrum (Fig. |2]i is well- 
described by a power law <p = <po (E / 1 TeV)~ r with a spectral 
photon index T = 2.0 + 0.1 sta t ± 0.2 sys and a flux normalization 
at 1 TeV of 0o = (4.2 + 0.8 stat ± 1.0 sys ) x 10~ 12 cirT 2 s _1 TeV" 1 . 
The integral flux F(l-10 TeV) = 3.8 x 10~ 12 ph cirT 2 s -1 is 
~17% of the Crab Nebula flux in the same energy range. The 
extracted flux points from the extended emission and the fitted 
power law are shown in Fig. [2] The results presented above have 
been cross-checked, using an independent calibration of the raw 
data and an alternative analysis chain. The cross-checks included 
a spectral analysis using the reflected region background method 
dBerge et al.ll2.007l) . which requires observations to be centered 
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outside of the emissive region and thus used only half of the 
available dataset. All cross-checks confirmed the primary results 
within the stated statistical uncertainties. 

The most recent observations and analysis by CANGAROO- 
III also give an indic ation of extended em ission in the vicin- 
ity of PSRB1706-44 dEnomoto et al.ll2009l) . However, their re- 
sults differ significantly from those given in this paper. For ex- 
ample, the morpholo gy of the VHE y-ray excess reported by 
lEnomoto et al.l d2009l) . using an ON-OFF background technique, 
is that of a source centered roughly at the pulsar position, as 
opposed to HESS J1708-443., whose centroid is clearly offset 
from the pulsar. Furthermore, CANGAROO-III measures a Crab 
Nebula-level integral flux (above 1 TeV) within 1.0° of the pul- 
sar, which is inconsistent with the ~18% Crab flux measured 
by H.E.S.S. in the same energy range. The difference is possi- 
bly due to the exact methods used for background subtraction; 
in the H.E.S.S. analysis, the OFF data are normalized to source- 
free regions of the ON data, because the background can vary 
significantly depending on the observing conditions. 

4. Origin of the VHE y-ray Emission 

While a superposition of a relic PWN created by PSR B 1706-44 
and SNR G 343. 1-2.3 cannot be excluded, each of these ob- 
jects individually could account for the observed VHE y-ray 
emission. The possible associations with HESS J1708-443. will 
be discussed in the following sections and both leptonic and 
hadronic scenarios will be considered. 



4.1. A Relic Nebula from PSR B1 706-44 

The pulsar PSR B 1706-44, which has a high spin-down lumi- 
nosity E = 3.4 x 10 36 erg s -1 , is energetic enough to power 
the observed VHE y-ray emission, which has a luminosity be- 
tween 1 and 10 TeV of L y « 9.9 x 10 33 (D/2.3kpc) 2 erg s~ l . 
The apparent conversion efficiency from rotational energy to y- 
rays in this energy range can be defined as ei-ioTeV = L y / E 
and for this case is ~0.3%, compatible with the efficien- 
cies (< 10%) of other VHE y-ray s ources which have well- 
established associations with PWN (!Gallantll2007l) . The pro- 
jected size of HESS J1708-443 corresponds to a physical size 
of ~12 (D/2.3kpc) pc (68% containment radius). These char- 
acteristics suggest a possible association between the VHE y- 



PWN/VHE associations, e.g. Vela X ( 


Aharonian et al.ll2006b) 


and HESS J 1825- 137 dAharonian et al. 


2006d). 



In a leptonic scenario, the VHE y-radiation originates from 
accelerated electrons which up-scatter ambient photons to VHE 
y-ray energies via inverse Compton (IC) scatteri ng. Compared 
to the size of the PWN in the radio (radius ~1.5') dGiacani et alJ 
12001 and the "bub ble" nebula seen in X-rays (radius ~1.8') 
dRomani et"ai1l2005h . the VHE y-ray PWN (sometimes referred 
to as a TeV PWN) would be a factor of ~10 larger. Similar 
differences in size have been obs erved in other TeV PWN as- 
sociations, e.g. HESSJ1825-137 dAharonian et al.ll2006dl) . and 
can be explained by the different energies, and hence cooling 
times, of the electrons which emit the X-rays and VHE y-rays. 
Assuming the magnetic field is uniform and that the average 
wind convection speeds in th e y-ray and X-ray emitti ng zones 
are both constant and similar, lAharonian et al.l (12005 dl) estimate 
the ratio of sizes 



^1-4 B 
Rx \ 



IOjuG 



Eh 



:eV 



where E^ e y is the mean energy in X-rays (2 keV) and £iev is the 
mean energy in VHE y-rays (0.9 TeV). However, in contrast to 
the PWN of PSR J 1826- 1334, where a magneti c field strength 
B - 1 uG was inferred from X-ray observations dGaensler et al.l 
l2003h . iRomani et al. I J2005I) estimated a magnetic field B as 



strong as 140^° uG within the 110" radius X-ray PWN of 
PSR B 1706-44, assuming the spectral break between the ex- 
trapolation of radio and X-ray spectra is due to radiative cool- 
ing of electrons. In such a high magnetic field, electrons that 
emit keV X-rays have comparable energies to those that emit 
TeV y-rays and therefore have comparable cooling times as well. 
Thus, the TeV PWN should be approximately the same size as 
the X-ray PWN, i.e. it should appear point-like considering the 
~5' H.E.S.S. PSF Furthermore, given that the ratio of X-ray to 
VHE y-ray energy flux is determined by the energy density in 
magnetic fields and IC target photon fields (only the cosmic mi- 
crowave background (CMB) is considered here), 



* 0.1(0.15-6)- 

^x 



(2) 



where B = 10~ 6 B_ 6 G dAharonian et al.l 1 19971) . the observed 
X-ray flux Fx = 3 3 x 10~ 13 erg crrT 2 s" 1 at £ k ev = 1.2 keV 
(lBeckeretal.1119951) can be used to predict the y-ray flux F y at 



> B -\ 



1 .4 TeV dde Jager & Diannati-Atai 



I2008h 



^TeV - 167£keVA _ ( 

assuming the value of B estimated bv lRomani et al.l d2005l) . This 
results in a predicted F y (lA TeV) = 1.7 x 10~ 16 erg crrT 2 s _1 , 
well below the level observable by H.E.S.S. Conversely, the ab- 
sence of VHE y-rays from the compact nebula (c.f. H.E.S.S. UL 
from Region A in Sect. 3) can be used together with Fx to calcu- 
late a lower limit on the magnetic field using Eq. 2. The resulting 
li mit, B ^ 2.5 uG, is con sistent with the magnetic field estimated 
bv lRomani et al] d2005l) . 

One way to reconcile the difference in emission region size 
and the high flux of the VHE y-ray emission is to assume that 
the size of the X-ray PWN is primarily governed by the extent 
of the high B-field region and that the magnetic field decreases 
by a large factor beyond the X-ray PWN. The electrons can then 
escape from the high B-field region and, by accumulating over 
a significant fraction of the lifetime of the pulsar, form a larger 
nebula which is visible only in VHE y-rays. The synchrotron 
cooling time of electrons that up-scatter CMB photons to ener- 
gies E y is given by 



^"synch 



40 



\140/iG/ I TeV/ 3 



(3) 



(1) 



In the 140 fiG field inside the X-ray PWN, the cooling time 
of up-scattered electrons producing 1 TeV y-rays is ~40 yr. 
Assuming a dominantly advective, rather than diffusive, trans- 
port process, the average flow speed needed to drive elec- 
trons from the pulsar position to the edge of the X-ray PWN 
(r ~ 110") within 40 yr is 0.1(D/2.3kpc) c. The im- 
plied flow speed i s reas onable following the arguments of 
iKennel & Coronitil d!984l) . although their model considers the 
case of the symmetric Crab Nebula, which is admittedly a sim- 
plification of the asymmetric PWN considered here. If the mag- 
netic field within the X-ray PWN was much higher than 140 /iG 
in the past, when most of the electrons were emitted, the restric- 
tions on the flow speed would become more stringent. However, 
in the low B-field region outside the X-ray PWN, the synchrotron 
lifetime increases. Even for a magnetic field strength of 10 pG, a 
value about three times as large as the interstellar magnetic field, 
the cooling time of the aforementioned electrons is about 8 000 
yr, almost half of the characteristic age of the pulsar (17 500 yr). 
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Fig. 3. Large field-of-view (FoV; 5.1° x 5.1°) VHE y-ray im- 
age of the region containing HESS J1708-443. The Gaussian- 
smoothed (<x = 0.10°) VHE y-ray excess from Fig. [TJ(left) is 
shown in color. The white contours indicate the intensity of the 
330 MHz radio emission detected with Very Large Array (VLA) 
observations (see also the green contours in Fig. [TJ left). The 
dotted white box represents the FoV covered by the VLA ob- 
servations. Outside of this region, green contour s indicate the 
lower -resolution 2.4 GHz radio continuum data dDuncan et al. 
1 1995b taken with the Parkes telescope; these observations have a 
half -power beamwidth of ~10.4'. This image also shows the ap- 
proximately ring-shaped region used for normalizing the back- 
ground in the ON-OFF background method (see Sect. 2.3); this 
region is delimited by the large dashed white circles, exclud- 
ing the known TeV source HESS J1702-420 located toward the 
Northwest. The Galactic plane is also located toward the NW 
and is indicated by a thick black dotted line. 



This relic TeV PWN scenario does not, however, explain the 
asymmetric morphology of HESS J1708-443, in particular its 
offset from the pulsar location, nor does it explain the lack of 
detectable VHE y-ray radiation from the location of the pul- 
sar itself, assuming that the pulsar and X-ray PWN are em- 
bedded in an extended shell of relic electrons. Such asymme- 
tries have been observed previously i n other TeV PWNs, e.g . 
HESSJ1718-385, HE SS J1809-193 (lAharonian et all l2007ah 
and HESS J 1825- 137 (lAharonian et al.ll2005dl l2006dl) . These 



asymmetries could be accounted for in two ways: as a direct re- 
sult of a high proper motion of the pulsar or as a result of a den- 
sity gradient in the ambient medium. The density gradient could 
lead to an asymmetry in the reverse shock of the supernova, or it 
could lead to a diff erent expansion velo c ity for the TeV y-ray 



In the preceeding discussion, it was assumed that the pul- 
sar dominantly accelerates electrons. If a consider able fraction 
of the a ccelerated particles a r e instead hadrons (e.g . iHorns et all 
(120061) : lAmato et all (120031) : iBednarek & Bartosikl (120031) ). the 
constraints imposed by the large magnetic field within the X-ray 
PWN are removed. In a hadronic scenario, n° mesons are pro- 
duced by inelastic interactions between accelerated protons and 
the ambient gas; they then decay, emitting VHE y-ray photons. 
In such a scenario, the VHE y-ray emission would trace the dis- 
tri bution of the target material. The bright radio arc, interpreted 
bv lBock & Gvaramadzel d2002l) as the compressed outer bound- 
ary of the former wind-blown bubble, could act as such a target 
due to its enhanced density, thereby also explaining the spatial 
coincidence with the H.E.S.S. source. Since the proton interac- 
tion time is long compared to the age of the pulsar, and assuming 
that the escape of protons from the region is sufficiently slow, all 
protons accelerated since the birth of the pulsar can contribute 
to the y-ray emission. However, to account for the high lumi- 
nosity of the VHE y-ray emission, the pulsar must have a high 
rotational energy and must efficiently convert rotational energy 
into proton acceleration. The total energy in accelerated protons 
W p in the energy range 10-100 TeV which is necessary to pro- 
duce the observed y-ray luminosity L y can be estimated from the 
relation 



W p (10 - 100 TeV) * T y x L y {\ - 10 TeV) . 



(4) 



where r y « 5 x 10 15 (n/cm 3 ) 1 s is the characteristic cool- 
ing time of protons through the n° production channel. The to- 
tal energy within the entire proton population Wp(tot) « 3 x 
10 49 erg(«/cnr 3 ) _1 (Z)/2.3kpc) 2 is then estimated by extrapo- 
lating the proton spectrum down to 1 GeV assuming the same 
spectral shape as the VHE y-ray spectrum, i.e. a power law with 
index T = 2.0. Assuming that a fraction rj of the pulsar's rota- 
tional energy E mt is converted into the energy within the proton 
population Wp(tot) = n E rot , then 



(cm- 3 )*°- 2 (2.3kpc) n '(lOms) ' 



(5) 



where E mt — and I « 1 x 10 45 g cm 2 is the moment of 

inertia of the pulsar. For a distance D = 2.3 kpc and an effi- 
ciency r] = 0.3, the initial rotation period Pq has to be as small 
as 6-12 ms for the ambient medium density to be in the range 
n * 1-5 cirT 3 . Although pulsars are thought to be born with ini- 
tial periods which are considerably shorter than their present pe- 
riods, the initial rotation period implied for PSRB 1706-44, in 
the above hadronic TeV PWN scenario, is even smaller than that 
of the Cr ab pulsar, the only case for w hich Pq is well-determined 
(19 ms) dManchester & Tavlorll 19771) . 

The hadronic PWN scenario is further disfavored by con- 
straints on the proton escape time. Under the common assump- 
tion that the proton diffusion coefficient is energy-dependent, i.e. 



emitti ng electrons (Bl ondin et al. 20011: van der Swaluw et all 



120011) . Simulations bv lvan der Swaluw et al.l(l200ll) demonstrate D(E p ) = D (E p / 10 GeV) 5 



that a displaced PWN can indeed be well-separated from its pul- 
sar. These explanations are in principle applicable to the case 
of HESS J1708-443; however, the pulsar's measured scintilla- 
tion velocity, less than 100 km s , renders the first explanation 
unlikely. The latter explanation favors a TeV PWN which is off- 
set toward a low density region. The available H I line emission 
data (see Fig. |4] (left) and the subsequent section) suggest that 
this might be the case, although it is not clear given the complex 
H I morphology. 



(6) 



with a power-law index 6 ~ 0.5, where Do is the diffusion coef- 
ficient at 10 GeV, one can estimate Do required to contain pro- 
tons with energy E p = 100 TeV within a certain distance of the 
pulsar after t = r c , since the diffusion radius R^f = 2 ^D{E)t 
for timescales less than the proton energy loss time, t «c r y . 
This containment region can be estimated to have an angu- 
lar size of ~0.7°, which is the approximate distance between 
the pulsar and the farthest significant VHE y-ray emission. At 
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Right Ascension (J2000) velocity [km/s] 



Fig. 4. Left: The grey-scale image shows the intensity of HI line emission in units of K k m s , measured by the Parkes radio 
telescope during the Southern Galactic Plane Survey (SGPS) (Mc Clure-Griffiths et al.ll2005h . The intensities are integrated in the 
velocity range -13.79 km s to -21.21 km s (shown as a shaded region in the velocity profile, right), corresponding to a near/far 
kinematic distance of 1.6-2.3 kpc / 13.0-13.7 kpc. Contours of the Gaussian-smoothed (cr = 0.10°) VHE y-ray excess are shown 
in blue. The red contours depict the intensity of the radio emission measured by the Very Large Array (VLA) at 330 MHz (see also 
Fig. Q] left). The radio data have been smoothed with cr = 0.03°. The white circle illustrates the integration region for the velocity 
profile shown on the right. Right: Velocity profile of H I line emission intensity, integrated over the region enclosed by the dashed 
circle on the left. The velocit y resolution is 0.0 8 km s _1 . The kinematic distance, shown in red, is derived from the velocity using 
the Galactic rotation curve o f lFichetal.l(fl989l 



the assumed pulsar distance, this region would have a physical 
size of 28 (D / 2.3 kpc) pc. The required diffusion coefficient 
Dq « 3.4xl0 25 {D 1 2.3 kpc) 2 cirrs -1 is found to be prohibitively 
low by a factor of 10-100 and can only be reconciled by assum- 
ing a very weak energy dependence, 5 < 0.2. 

4.2. SNR G 343.1 -2.3 

The VHE y-ray source HESS J 1708 -443 is partially coincident 
with the bright radio arc and the surrounding diffuse emission 
of the SNR, visible in the 330 MHz observations taken with the 
VLA (see contours in Fig. Q] left, and Fig. IA.U . The centroid 
of the H.E.S.S. source is consistent with the apparent center of 
the brig ht radio arc (g?.on o = 17 h 8 m , <5 2 ooo = -44°16'48"; as de- 
fined in lAharonian et all (12005 al) . The extension of the VHE y- 
ray excess (68% containment radius: 0.29° ± 0.04°) i s compati- 
ble wi th the the radius of the radio shell (~0.27°) fit bv lFrail et al.l 
(fl99l using VLA observations at 90 cm. The 95% contain- 
ment radius (0.71°; used for spectral extraction) of the y-ray 
excess completely encloses the radio shell, whos e approximate 
bound ary was estimated at a radius of ~0.42° by iRomani et alj 
(120051) using ATCA (Aus tralia Telescope Compa ct Array) ob- 
servations at 1384 MHz dDodson & Golapll2002l) . Thus, while 
the majority of the VHE y-ray emission is located within the 
radio shell, emission from the shell itself cannot be excluded. 
Due to low statistics in the current VHE dataset, no further 
conclusions can be made regarding morphological similarities. 
No significant VHE emission was detected from the spatially- 
extended, diffuse emission visible farther to the Southeast of the 
bright radi o arc, seen in the low -resolution 2.4 GHz continuum 
radio data dDuncan et al ] 11995b shown in Fig. [3] although the 
offset-corrected exposure in this region is very low (between 
~4 and 10 h) since all the H.E.S.S. observations were centered 
near the pulsar. This diffuse radio emission was interpreted by 



iBock & Gvaramadzel d2002l) as originating from the eastern half 
of the expanding SNR shell, propagating into a low-density re- 
gion. 

Similar to the potential association with the PWN of 
PSRB 1706-44, both leptonic and hadronic scenarios will be 
considered for VHE y-ray production. The leptonic scenario suf- 
fers from the non-detection of the SNR at X-ray energies. The 
VHE y-ray spectrum is hard and extends up to 20 TeV; assum- 
ing IC scattering in the Thomson regime, the electrons which 
up-scatter CMB photons to 20 TeV have an energy of ~80 TeV. 
For a reasonable magnetic field strength of 5 jt/G, these elec- 
trons would emit synchrotron photons with an energy of ~ 1 ke V, 
i.e. photons within the detectable energy range of current X- 
ray telescopes. Unfortunately, this p rediction cannot be tested 
because the X-ray UL calculated by iBecker et al.l (1 19951) using 
ROSAT was derived from a relatively small part of the shell; no 
stringent UL on the X-ray flux from an extended region within 
0.7° of the H.E.S.S. source can be derived (W. Becker, per- 
sonal communication) due to its large extension and the vicinity 
of the luminous low -mass X-ray binary (LMXB) 4U 1705-440 
dBecker et al.ll 1995b . whose stray light may be obscuring diffuse 
X-ray emission from the SNR. It is also possible that the X-ray 
emission is inherently weak and can not be detected due to the 
relatively high interstellar absorption (IBecker et al.l 1995b . 

In the hadronic scenario, synchrotron radiation is expected 
only from secondary electrons, and the lack of X-ray detection 
can easily be accounted for. Assuming a total energy of 10 51 erg 
is released in the supernova explosion, an acceleration efficiency 
of e = 0.15 and a distance D = 2.3 kpc, an average proton den- 
sity of n « 1.5 cirT 3 - a value slightly larger than the average 
Galactic ambient density - is sufficient to explain the previously 
estimated (Sect. 4.1) energy content of Wp(tot) w 3 x 10 49 erg 
within the proton population. 
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Given the various scenarios that have been proposed to ex- 
plain the origin of the bright radio arc, there are many differ- 
ent possibilities as to how the SNR could be associated with 
HESS J1708-443. In one scenario, the SNR G 343.1-2.3 is ex- 
panding symmetrically into the interstellar medium (ISM), and 
the intensity variations which form the radio arc are due to local 
density differences in the ISM. An association between the SNR 
and the pulsar PSRB 1706-44, a controv ersial scenario which 
is stil l debated in the comm unity (see e.g. iBock & Gvaramadzd 
l2002t iRomani et all l2005h . would make the SNR rather old, 
(9(10000 yr), and place it in the late Sedov-Taylor phase, or, 
more likely, in the radiative phase. If the SNR is in the radia- 
tive phase, the ambient material swept up by the SNR should be 
visible in CO or H I data. Unfortunately, no high-resolution CO 
data are publicly available at the moment, but there is evidence 
for a ring-like structure in the H I line emission survey of the 
Parkes telescope, as shown in Fig. |4] (left). The structure is best 
visible in the velocity range -13.79 km s to -21.21 km s _1 
corresponding to a near/far kinematic distance of 1.6-2.3 kpc 
/ 13.0-13.7 kpc (Fig. [4] right). The near distance is compati- 
ble with the pulsar distance (2.3 kpc) A rough estimate of the 
mass of the H I structure, extracted from the circular region in 
Fig. |4] (left), is ~6 x 10 3 Mq. Assuming radial symmetry, this 
corresponds to an original density of the swept-up mass of a few 
protons cirT 3 , comparable to the density requirements imposed 
by the observed y-ray flux. The bright radio arc and the VHE 
y-ray emission spatially coincide with only one half of the H I 
shell-like structure. This morphology could arise because of the 
additional dependence of the radio and y-ray emission on the 
target density, which is likely larger closer to the parent MC. 

The HI shell has a radius of ~0.4°, which, assuming a 
distance of D = 2.3 kpc, corresponds to a physical radi us of 
~16pc. Following t he app roaches of ICioffi et al] ( 1 1988b and 
iTruelove & McKed d 19991) and further assuming an age of 
17 000 yr, a 10 Mq progenitor star, and an energy release of 
10 51 erg, the ambient density necessary to explain the size of 
the HI shell is n ~ 0.7 cm ' 3 and the resulting shock vel ocity 
is -400 km s" 1 . Following fPtuskin & Zirakashvilil d2005l) . the 
maximum proton energy attainable is then (9(10 TeV), likely too 
low to explain the observed TeV emission, which extends up to 
20 TeV. The spectral energy distribution (SED) of y-rays pro- 
duced in the interactions of mono-energetic protons (and subse- 
quent decay of pions) drops sha rply beyond roughly 15% of the 
original proton energy, see e.g. dKelner et al.ll2006l) . Therefore, 
the parent proton population giving rise to the observed VHE y- 
ray emission should extend up to about 100 TeV, a limit which 
is - as the example calculation above illustrates - increasingly 
difficult to explain as the age of the system increases. Indeed, 
the y-ray emitting SNR shells which have bee n unambiguously 
identifi ed so far, such as RX J1 713.7-3946 ( Aharonian et al.l 
I2007bl) and RX J0852.0-4622 dAharonian et al.l l2007d) . "are 
much younger (~2 000 yr). 

The aforementioned constraints are removed if the SNR ex- 
pands first into a bubble blown by the progenitor star's wind into 
the ISM. Due to the low density inside the wind-blown bubble, 
the velocity of the expanding shock is much higher than an- 
ticipated and protons can be accelerated to very high energies. 
When the shockfront reaches the outer boundary of the wind- 
blown bubble, the high-energy protons are released to interact 
with the dense environment outside of the bubble and produce 
VHE y-rays in the p r ocess. In this scenario, first proposed by 
Bock & Gvaramadzd d2002l) . the bright radio arc is created by 
the former boundary of the wind-blown bubble which has been 
overtaken by the expanding SNR shockfront. The offset of the 



pulsar position from the center of the radio arc does not hin- 
der the association between the SNR and the pulsar, since the 
progenitor star, whose wind has produced the bubble, can have 
traversed the bubble's boundary during its evolution, before it 
became a supernova. The constraints on the VHE y-ray produc- 
tion imposed by the large implied age of 17 000 yr do not apply 
in this case since the protons now interacting within the dense 
ambient medium to produce y-rays could have been acceler- 
ated in the past, when the SNR shock velocity was still high. 
However, this would require an extremely low diffusion coeffi- 
cient (Do ~ 2 x 10 25 cm 2 s ), similar to the case of the 
hadronic PWN scenario (see Sect. 4.1). 

This discussion of a putative association between 
HESSJ1708-443 and the SNR G 343.1-2.3 is based on 
the assumption that the SNR and the pulsar PSRB 1706-44 
were created at the same time. If this assumption proves 
to be wrong, then very little is known about the SNR. The 
age estimate using a Sedov-Taylor model is about 5 000 yr 
dMcAdam et al.l 119931: iNicastro et al.l 119961) . The younger age 
would further ease the proton acceleration to energies beyond 
100 TeV. 

To summarize, the radio emission from SNR G 343. 1-2.3, 
which may originate from the interaction of the SNR with an 
ambient MC, is partially coincident with HESS J1708-443, sug- 
gesting a plausible association which could account for at least 
part of the VHE y-ray emission observed. However, the puta- 
tive associations between the SNR and the pulsar or between 
the SNR and the shell-like structure discovered in HI suggest 
that the SNR is in a later evolutionary stage than other SNRs 
previously-detected in the VHE regime. 

4.3. Other Nearby Celestial Objects 

There are other celestial objects nearby, i.e. within the emission 
region of HESS J 170 8-443, notably the LMXB 4U 1705-440 
dForman et al l 1 19781) and the radio source PMN J 1708 -44 19 
i Wright etalT 1 19941) The LMXB is a well-studied type 1 



burster dSztaino et a l.ll985l) located at g 7n m = 17 h 08 m 54.46 s and 
<5 2 ooo = -44°6'7.35" dDi Salvo et1ii1l2005l) . i.e. it is offset from 
the centroid of the VHE emission by 0.25°. Considering this off- 
set and the extended nature of the VHE y-ray source, an associa- 
tion is highly unlikely since an X-ray binary would appear point- 
like to H.E.S.S. Theoretical predictions for VHE y-ray emis- 
sion from LMXBs focus on those with relativistic jets (micro- 
quasars); 4U 1705-440 does not exhibit jets. Furthermore, no 
LMXBs have been detected in the VHE y-ray regime, despite 
the extensive coverage of the H.E.S.S. Galactic Plane Survey. 

The radio source PMN J 1708 -44 19 is located at 
g 7 ooo = 17 h 08 m 30 s and <5 2 ooo = -44°19'07" dWright et al.l 
Il994l) . The local maximum in the radio contours at the center 
of Fig. [T] is largely due to this very bright point source, clearly 
visible in the 330 MHz VLA radio image (see Fig. IA. lb . 
Although its position is compatible with the centroid of the 
H.E.S.S. source, an association between the two is unlikely 
given the spectrum of the radio source. Using data from the 
VLA (at 330 MHz and 1.4 GHz), Molonglo Galactic Plane 
Survey (MGPS; cataloged as J 170828-441823 at 840 MHz, 
and Parkes-MIT-NRAO (PMN; cataloged as PMN J 1708-44 19 
at 4.8 GHz), we derive a spectral index a = -0.81 ± 0.08 in 
the radio dom ai n, con sistent with the value a = -0.9 derived 
by iFrail et all d 19941) over a narrower range in frequency, 
from 330 to 840 MHz. The steep spectral index suggests that 
PMN J 1708-44 19 is extragalactic, since Galactic point-like 
sources are typically compact HII regions, for which the 
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radio spectral index is positive; therefore, it is unlikely to be 
associated with the extended emission of HESS J 1708-443. 
Upon a deeper inspection of this source using high-resolution 
unpublished ATCA radio data (Dodson & Golap, personal 
communication), this bright source can be further resolved 
into two sources. However, the spectral indices above were 
calculated on the basis of observing it as a single unresolved 
source, because this is the way that the low-resolution radio 
surveys detected them. 

In order to quantify the contribution any putative, unre- 
solved point source could make to the flux observed from 
HESS J1708-443, one can compare the symmetric 2D Gaussian 
curve of a point source to that of the extended H.E.S.S. source. 
This demonstrates that any single unresolved point source 
could not account for more than ~6% of the total flux from 
HESSJ1708-443. 



5. Summary 

H.E.S.S. observations of the y-ray pulsar PSRB 1706-44 have 
led to the detection of an extended (<x = 0.29° + 0.04°) source 
of VHE y-ray emission, HESS J1708-443, in the Galactic plane. 
Its energy spectrum is well-described by a power law with a pho- 
ton index F = 2.0 + 0.1 stat ± 0.2 sys and a normalization at 1 TeV 
of 0o = (4.2 ± 0.8 sta , + 1.0 sys ) xl0~ 12 cm" 2 s" 1 TeV 1 . The cor- 
responding integral flux F( 1-10 TeV) = 3.8 x 10 _12 phcm 2 s _1 
is roughly 17% of the Crab Nebula. The possible associations of 
HESS J1708-443 with an offset, relic PWN of PSRB1706-44 
and with the partial shell-type SNR candidate G 343.1-2.3 have 
been discussed using additional radio and H I line emission data. 
Given the extended nature of the TeV source and the limited 
statistics, it is unclear if the emission is associated with the 
PWN, located at the edge of the H.E.S.S. source, or with the 
SNR, in which the pulsar is thought to be embedded. Based 
on energetics and a wealth of information at other wavelengths, 
neither interpretation can be excluded at this time; furthermore, 
the possibility remains that both sources contribute to the to- 
tal observed VHE y-ray emission. High-spatial-resolution CO 
mapping of this region would improve our understanding of the 
molecular environment and might help to identify a preferred 
MWL counterpart to HESS J 1708-443. Deeper exposure in the 
TeV regime would also provide vital statistics and enable more 
detailed morphological and spectral studies. 
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Appendix A: Comparison with H.E.S.S. 2003 
dataset 

A. 1. Recalculation of upper limits using the 2003 dataset 

In its commissioning phase, the H.E.S.S. IACT observed the 
region around th e energetic PSRB1706 -44 between April 
and July 2003 (lAharonian et alJ l2005al) . No evidence for 
statistically-significant VHE y-ray emission was found at the 
pulsar position nor from a region encompassing the partial shell- 
type SNR G 343. 1- 2.3. Upper limits (ULs) to the integral flux 
were published in (lAharonian et al 1 l2005ah . The integral flux 
now measured by H.E.S.S. (see Sect. 3) is not compatible with 
those originally-published ULs, a discrepancy which motivated 
a re-analysis of the 2003 H.E.S.S. dataset for this region, using 
the current H.E.S.S. software. 

Although H.E.S.S. is currently an array of four IACTs, it was 
operating as a two-telescope array from February to December 
2003. The 2003 observations yielded a dataset with a livetime of 
14.3 h, an average zenith angle of ~26°, and an energy thresh- 
old, estimated from Monte Carlo simulations, of ~350 GeV. For 
y-hadron separation, standard cuts were used, which require a 
minimum of 80 p.e. to be recorded per shower image. 

Integral flux ULs were calculated from three a priori de- 
fined circular regions: a Standard point-like (8 = 0.14°) re- 
gion centered at the position of PSR B 1706-44, a CANGAROO- 
like region (8 = 0.22°) also centered at the pulsar position, 
and a Radio arc region centered at the apparent center of the 
SNR G 343.1-2.3 (a 7 tm = 17 h 08 m , <5 2 ooo = -44°16'48"; as de- 
fined in lAharonian et al ] (l2005ah V with a radius 8 = 0.60° in 
order to completely enclose the complex radio structure. The 
CANGAROO region is disregarded for the remainder of this ap- 
pendix, because its sole purpose was to compare the H.E.S.S. re- 
sult with the ori ginal CANGAROO-I de tection, which has since 
been rescinded (lYoshikoshi et al.l l2()09i. and focus primarily on 
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Fig. A.l. Integral flux upper limits (ULs) from the Radio arc re- 
gion encompassing both PSRB 1706-44 and SNR G 343.1-2.3. 
The red po ints represent the previou sly-published integral flux 
ULs from lAharonian etail (12005 al) using the 2003 H.E.S.S. 
dataset. The black solid line is based on a revised calculation 
of the UL using the same dataset. 
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Table A.l. Event statistics and background parameters for the 
analyses of the Radio arc region around PSRB 1706-44 and 
SNR G 343. 1-2.3. Row 1 gives the integral flux upper lim- 
its (99% confidence level) from both analyses. The number of 
events N in the circular (radius 9 = 0.6°) on-source (ON) and 
ring-shaped off-source (OFF) regions are given in rows 2 and 
3, the normalization factor a (the ratio of O N to OFF area) i n 
row 4, excesses and significances (according to lLi & Mai (1 19831) 1 
in rows 5 and 6, and the ring parameters in rows 7 and 8. The 
statistics and upper limits presented here were obtained using 
the 2003 H.E.S.S. dataset only. 



the Radio arc region (equivalent to Region B; see Sect. 3), which 
is very similar to the region from which extended VHE y-rays 
are now detected (Region C; see Sect. 3). 

Background subt raction was perf ormed using the Ring 
Background Method dBerge et alj|2007l) in both the original and 
revised analysis. The exact inner and outer ring radii, ri nner and 
/"outer respectively, used in the original analysis were not docu- 
mented; however, the inner ring radius is typically chosen to be 
slightly larger than the on-source (ON) region (radius 6 = 0.60°) 
and the normalization factor a = 1/7 (the ratio of the ON to 
off-source (OFF) area). Therefore, for the re-analysis, the in- 
ner ring radius is chosen to be 0.65°, which, given a, leads to 
'"outer = 1 -35°. The number of events in the ON and OFF regions, 
AfpN and A^off, respecti vely, is found to match those given in 
lAharonian et al.l (12005a) to within 8% (see Table lATTT i. demon- 
strating that the ring parameters adopted in the re-analysis are 
approximately equal to those in the original analysis. No ex- 
clusion region was placed on the now known-to-exist source, 
HESS J1708-443, i.e. the source is not excluded from OFF re- 
gions. In practice this has a negligible effect, since there is lit- 
tle y-ray emission from HESS J1708-443 beyond r; nner = 0.65° 
from its centroid. 

Th e UL (99% confidence level using iFeldman & Cousinsl 
(fl998h on the integral flux from the Radio arc region (Region 
B) was originally found to be F{> 0.35 TeV) < 5.8 x 
10~ 12 ph ctrT 2 s , equivalent to ~5% Crab, assuming the spec- 
trum is describe d by a power law with a spectral index T = 
2.5 {Method A in lAharonian et al.l (1200531) 1. An alternative UL, 
F{> 0.50 TeV) < 3.5 x 10" 12 ph ctrT 2 s~\ also equivalent to 
~5% Crab, was calculated us i ng a m ethod {Method B, described 
in detail in lAharonian et al.l (12005 al) l which made no assump- 
tions concerning the source spectrum. These ULs are shown in 
Fig. IA.ll where they are compared to the revised calculation (us- 
ing Method A) of the integral flux UL, plotted as a function of 
threshold energy E. The revise d UL is clearly h i gher (le ss strin- 
gent) than the one published in lAharonian et al.l d2005al) . For ex- 
ample, the integral flux above 0.35 TeV is F{> 0.35 TeV) < 9.7 x 
10~ 12 ph cirT 2 s , equivalent to 9% Crab and above 0.50 TeV 
is F{> 0.50 TeV) < 7.6 x 10~ 12 ph cirT 2 s" 1 , equivalent to 
12% Crab, again assuming F = 2.5. See Table |A~T1 for a sum- 



mary and comparison of the event statistics and other analysis 
parameters from both analyses. 

The use of two-telescope data resulted in a lower sensi- 
tivity at the time but would not have had any negative im- 
pact on the original determination of ULs from the vicinity 
of PSRB 1706-44. H.E.S.S. currently uses a stereo trigger im- 
plemented at the hardware level to select extended air show- 
ers (EASs) simul taneously detected by at least two telescopes 
dFunk et al.ll2004l) . However, from February to July 2003, when 
the original observations of PSRB 1706-44 were carried out, 
it used an off-line triggering mode, since the central hardware 
trigger had not yet been installed. In software stereo mode, 
each recorded EAS receives a time stamp via a GPS (Global 
Positioning System) clock. The time stamps are then used in the 
offline data analysis to identify EASs which were observed in 
coincidence by the two telescopes. The use of a software stereo 
trigger, while not as efficient as the hardware stereo trigger cur- 
rently in use, is not expected to have contributed significantly to 
the discrepancy between the original ULs and the new results. 

After investigating various possible reasons for the discrep- 
ancy, it remains unknown why the previously-determined ULs 
were so low, leaving human error or undocumented changes 
in the analysis software used at the time as possible explana- 
tions. It is important to note that many other published results 
based on data taken during H.E.S.S.'s commissioning phase have 
been subsequently confirmed by the full four-telescope array 
with a hardware trigger , e.g. obser vations of RXJ1713-3 946 
(lAharonian et al.ll2007bl) and Sgr A* (lAharonian et alj|2009l) . 

A.2. Compatibility between detected flux and 2003 upper 
limits 

The upper limits calculated in the previous section cannot be di- 
rectly compared to the new H.E.S.S. results, based on the 2007 
dataset (presented in Sect. 3), because they assume a spectral 
index F = 2.5 and a low energy threshold. The new VHE y- 
ray source, HESS J 1708-443, has a much harder spectral index 
F = 2.0 ± 0.1 s tat ± 0.2 sys . Furthermore, the minimum energy 
threshold of the H.E.S.S. array has increased due to the reflectiv- 
ity of the IACT mirrors diminishing from 2003 to 2007, which 
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reduces the array's ability to detect faint EASs initiated by lower- 
energy y-rays. 

A re-analysis of the 2003 dataset, using the current H.E.S.S. 
software and assuming F = 2.0, yields a flux upper limit (99% 
confidence level) F(> 0.6 TeV) < 6.3 x 10~ 12 ph cm" 2 s _1 , 
equivalent to -13% Crab, for Region B. Analysis of the 2007 
dataset shows a statistically-significant signal from Region B, 
F(> 0.6 TeV) ~ 6.5 x 10~ 12 ph cm" 2 s _1 , equivalent to -13% 
Crab. These two flux values are statistically compatible, given 
the typical uncertainties in the measured flux normalization 
(+ -20%) and spectral index (+ -0.2). 

Furthermore, the new analysis results for HESS J1708-443 
(Region C), based on the 2007 dataset, have been confirmed us- 
ing an independent data calibration and analysis chain, and the 
cross-check analysis is also compatible with the presented re- 
sults. 
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Right Ascension (J2000) 

Fig. A.l. Image showing the intensity of the radio emission measured by the Very Large Array (VLA) (Frail et a l.ll994l) at 330 MHz 
in the vicinity of PSR B 1706-44, smoothed with a Gaussian of width <x = 0.03°. The observations have a half-power beamwidth of 
0.03° x 0.015°. The radio arc of the partial shell-type SNR G 343.1-2.3 is clearly visible as well as the diffuse emission both inside 
and outside the arc. The white contours correspond to a smoothed VHE excess of 0.14, 0.17, and 0.21 y-rays arcmkT 2 , taken from 
the image in Fig.Q] left. The horizontal stripes visible at Dec = -44°40 / and Dec = -43°47.5' are imperfections which resulted 
from the joining of data to form the final wide-field image dFrail et al.lll994h . The bright point source at the center of the radio image 
is PMN J 1708 -44 19, likely an extragalactic object seen in projection (see Sect. 4.3). 



